Evolution of Second-Order Cosmological Perturbations and Non-Gaussianity 
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We present a second-order gauge-invariant formalism to study the evolution of curvature perturba- 
tions in a Friedmann- Robertson-Walker universe filled by multiple interacting fluids. We apply such 
a general formalism to describe the evolution of the second-order curvature perturbations in the 
standard one-single field inflation, in the curvaton and in the inhomogeneous reheating scenarios for 
the generation of the cosmological perturbations. Moreover, we provide the exact expression for the 
second-order temperature anisotropies on large scales, including second-order gravitational effects 
and extend the well-known formula for the Sachs- Wolfe effect at linear order. Our findings clarify 
what is the exact non-linearity parameter /nl entering in the determination of higher-order statis- 
tics such as the bispectrum of Cosmic Microwave Background temperature anisotropies. Finally, we 
compute the level of non-Gaussianity in each scenario for the creation of cosmological perturbations. 

PACS numbers: 98.80.cq; DFPD-A-03-36 



I. INTRODUCTION 

Inflation [1,2] has become the dominant paradigm for 
understanding the initial conditions for structure for- 
mation and for Cosmic Microwave Background (CMB) 
anisotropy. In the inflationary picture, primordial den- 
sity and gravity-wave fluctuations are created from quan- 
tum fluctuations "redshifted" out of the horizon during 
an early period of superluminal expansion of the uni- 
verse, where they are "frozen" [3-7] . Perturbations at the 
surface of last scattering are observable as temperature 
anisotropy in the CMB, which was first detected by the 
COsmic Background Explorer (COBE) satellite [8-10]. 
The last and most impressive confirmation of the infla- 
tionary paradigm has been recently provided by the data 
of the Wilkinson Microwave Anistropy Probe ( WMAP) 
mission which has marked the beginning of the preci- 
sion era of the CMB measurements in space [11]. The 
WMAP collaboration has produced a full-sky map of the 
angular variations of the CMB, with unprecedented ac- 
curacy. WMAP data confirm the inflationary mechanism 
as responsible for the generation of curvature (adiabatic) 
super-horizon fluctuations. 

Despite the simplicity of the inflationary paradigm, 
the mechanism by which cosmological adiabatic pertur- 
bations are generated is not established. In the standard 
scenario associated to one-single field models of inflation, 
the observed density perturbations are due to fluctua- 
tions of the inflaton field itself. When inflation ends, the 
inflaton oscillates about the minimum of its potential and 
decays, thereby reheating the universe. As a result of the 
fluctuations each region of the universe goes through the 
same history but at slightly different times. The final 
temperature anisotropies are caused by the fact that in- 
flation lasts different amounts of time in different regions 
of the universe leading to adiabatic perturbations. Un- 
der this hypothesis, the WMAP dataset already allows to 



extract the parameters relevant for distinguishing among 
single- field inflation models [12]. 

An alternative to the standard scenario is represented 
by the curvaton mechanism [13-16] where the final curva- 
ture perturbations are produced from an initial isocurva- 
ture perturbation associated to the quantum fluctuations 
of a light scalar field (other than the inflaton), the curva- 
ton, whose energy density is negligible during inflation. 
The curvaton isocurvature perturbations are transformed 
into adiabatic ones when the curvaton decays into radi- 
ation much after the end of inflation. Recently, another 
mechanism for the generation of cosmological perturba- 
tions has been proposed [17-19] dubbed the inhomoge- 
neous reheating scenario (sometimes called the modu- 
lated reheating scenario)*. It acts during the reheating 
stage after inflation if super-horizon spatial fluctuations 
in the decay rate of the inflaton field are induced during 
inflation, causing adiabatic perturbations in the final re- 
heating temperature in different regions of the universe. 

Contrary to the standard picture, the curvaton and 
the inhomogeneous reheating scenario mechanism exploit 
the fact that the total curvature perturbation (on uni- 
form density hypersurfaces) £ can change on arbitrarily 
large scales due to a non-adiabatic pressure perturbation 
which may be present in a multi-fluid system [21-25]. 
While the entropy perturbations evolve independently of 
the curvature perturbation on large scales, the evolution 
of the large-scale curvature is sourced by entropy pertur- 
bations. 



"The idea that the total curvature perturbation may be af- 
fected on large scales by entropy perturbations when there 
exists a scalar field affecting particle masses or couplings con- 
stants controlling the reheating process has first been sug- 
gested in [20]. 
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Fortunately, the standard and the curvaton scenarios 
have different observational signatures. The curvaton 
and the inhomogeneous reheating scenario allow to gen- 
erate the observed level of density perturbations with a 
much lower scale of inflation and thus generically predict 
a smaller level of gravitational waves. 

More interestingly, the various scenarios for the gener- 
ation of the cosmological perturbations predict different 
levels of non-Gaussianity which is usually parametrized 
by a dimensionless non-linear parameter /nl- The non- 
Gaussianity can be large enough to be detectable either 
by present CMB experiments, the current WMAP [26] 
bound on non-Gaussianity being |/nl| < 10 2 , or by fu- 
ture planned satellites, such as Planck, which have enough 
resolution to detect non-Gaussianity of CMB anisotropy 
data with high precision [27] . 

Since a positive future detection of non-linearity in 
the CMB anisotropy pattern might allow to discriminate 
among the mechanisms by which cosmological adiabatic 
perturbations are generated, it is clear that the precise 
determination of the non-Gaussianity predicted by the 
various mechanisms is of primary interest. 

In this paper we achieve different goals: 

• We present a gauge-invariant formalism at second- 
order to study the evolution of curvature pertur- 
bations in a Friedmann- Robertson- Walker universe 
filled by multiple interacting fluids. We apply this 
general formalism to describe the evolution of the 
second-order gauge invariant curvature perturba- 
tions from inflation down to the radiation/matter 
phase in the various scenarios for the generation of 
the cosmological perturbations. In particular, we 
study the evolution of the total curvature pertur- 
bation during the non-adiabatic phase when either 
the inflaton field (for the standard and the inhomo- 
geneous reheating scenario) or the curvaton field os- 
cillate and decay into radiation. We show that the 
total curvature perturbation is exactly conserved in 
the standard scenario during the reheating phase 
after inflation and is sourced by a non-adiabatic 
pressure in the curvaton and the inhomogeneous 
scenarios. 

• We present the exact expression for the second- 
order temperature fluctuations on super-horizon 
scales, extending the well-known expression for the 
Sachs- Wolfe effect at linear order. 

• We provide the exact definition for the non-linear 
parameter /nl entering in the determination of 
higher-order statistics as the bispectrum of the tem- 
perature anisotropies. 

• We compute the exact expressions for the non- 
linear parameter /nl in the different scenarios for 
the creation of cosmological perturbations. 

The paper is organized as follows. In Section II we study 
the time-evolution for the system composed by a generic 



scalar field decaying into radiation which can be used to 
study a reheating phase. In Section III we describe in 
a gauge-invariant manner the evolution of cosmological 
perturbations at linear order while in Section IV we ex- 
tend this analysis to second-order providing the evolution 
equations for the curvature perturbations. In Section V 
we apply our findings to the various scenarios for the gen- 
eration of the cosmological perturbations and give the ex- 
act expressions for the second-order total curvature per- 
turbation. In Section VI we compute the expression for 
the second-order temperature fluctuations on large scales 
and the corresponding non-linear parameter. In Section 
VII the predictions of the non-linear parameter for the 
various scenarios are provided. Finally, our conclusions 
are contained in Section VIII. 



II. THE METRIC AND THE BACKGROUND 
EQUATIONS 

In this section we provide the equations describing 
the time-evolution for the system composed by a generic 
scalar field and radiation. In particular, we will focus 
on the evolution of cosmological perturbations on large 
scales, up to second-order, for a system composed by a 
scalar field oscillating around the minimum of its po- 
tential and a radiation fluid having in mind the physi- 
cal case of reheating. Averaged over several oscillations 
the effective equation of state of the scalar field tp is 
w v = (P<p/p v ) = 0, where P v and p v are the scalar 
field pressure and energy density respectively. The scalar 
field is thus equivalent to a fluid of non-relativistic par- 
ticles [28]. Moreover it is supposed to decay into radia- 
tion (light particles) with a decay rate T. We can thus 
describe this system as a pressureless and a radiation 
fluid which interact through an energy transfer triggered 
by the decay rate T. We follow the gauge-invariant ap- 
proach developed in Ref. [25] to study cosmological per- 
turbations at first-order for the general case of an arbi- 
trary number of interacting fluids and we shall extend 
the analysis to second-order in the perturbations. In- 
deed the system under study encompasses the dynamics 
of the three main mechanisms for the generation of the 
primordial cosmological density perturbations on large 
scales, namely the standard scenario of single field in- 
flation [1,2], the curvaton scenario [13-16], and the re- 
cently introduced scenario of "inhomogeneous reheating" 
[17-19,29]. In particular our second-order results allow 
us to calculate the level of non-Gaussianity in the gravita- 
tional potential in the latter scenario, where the inflaton 
decay rate has spatial fluctuations <5r, produced when 
the couplings of the inflaton field to normal matter de- 
pends on the vacuum expectations value of some other 
light scalar fields x present during inflation. 

Let us consider the system composed by the oscillating 
scalar field ip and the radiation fluid. Each component 
has energy- momentum tensor T^K and T^K. The total 
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1 (<p) 



27^ is covariantly con- 



energy momentum T' 1 " 
served, but allowing for an interaction between the two 
fluids [30] 
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where Q( v ) an d are the generic energy-momentum 
transfer to the scalar field and radiation sector respec- 
tively and are subject to the constraint 
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The energy-momentum transfer Qfo and Q 
decomposed for convenience as [30] 

Q(<p) = Q<p uV + f(tp) > 



(7) 



(2) 
can be 
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where the /"'s are required to be orthogonal to the the 
total velocity of the fluid u v '. The energy continuity equa- 
tions for the scalar field and radiation can be obtained 



from UvV^T^ = u v Q v u 
hence from Eq. (3) 
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In the case of an oscillating scalar field decaying into ra- 
diation the energy transfer coefficient Q v is given by [20] 



(5) 



where T is the decay rate of the scalar field into radiation. 

A consistent study of the density perturbations re- 
quires to consider the metric perturbations around a 
spatially-flat Friedmann-Robertson- Walker (FRW) back- 
ground as well. In the following we have used the 
expression for the metric perturbations contained in 
Refs. [31,32]. The components of a perturbed spatially 
flat Robertson- Walker metric can be written as 
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where a(r) is the scale factor and r = J dt/a is the con- 
formal time. The standard splitting of the perturbations 
into scalar, transverse (i.e divergence- free) vector parts, 
and transverse trace-free tensor parts with respect to the 
3-dimensional space with metric Sij can be performed in 
the following way: 



(7) 



1&=Di i x {T) +diX { p+d j ^ + ^, (8) 

where r = 1, 2 stand for the order of the perturbations, u>i 
and Xi are transverse vectors (d l (j^ — d z x^ = 0), x\^ 
is a symmetric transverse and trace-free tensor (d l x\^ = 

0: X^ r ' — 0) and = didj — (1/3) Sij d k dk is a trace- 
free operator. Here and in the following latin indices are 
raised and lowered using 6 % i and Sij, respectively. 
For our purposes the metric in Eq. (6) can be simpli- 
fied. In fact, first-order vector perturbations are zero; 
moreover, the tensor part gives a negligible contribution 
to second-order perturbations. Thus, in the following 
we can neglect uj^\ x'P and X$ • However the same 
is not true for the second-order perturbations. In the 
second-order theory the second-order vector and tensor 
contributions can be generated by the first-order scalar 
perturbations even if they are initially zero [31]. Thus 
we have to take them into account and we shall use the 
metric 
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= -a 2 (r)(l + 20 (1) +4» (2) ) , 
= a 2 (r) [d^) + 1 - di ^) + 1 - 
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The controvariant metric tensor is obtained by requiring 
(up to second-order) that g^ v g vX = S^ A and it is given 

by 
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_i(aV (2) + 5V (2) + x ij(2) ) 

- d i wW& i uW - 4V> (1) D ij x (1) +D lk X (1) D\x (1) \ • (10) 



A. Background equations 

The evolution of the FRW background universe is gov- 
erned by the Friedmann constraint 
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and the continuity equation 

p' = -m ( P + p) 



(ii) 



(12) 



where a prime denotes differentiation with respect to con- 
formal time, H = a' /a is the Hubble parameter in con- 
formal time, and p and P are the total energy density 
and the total pressure of the system. The total energy 
density and the total pressure are related to the energy 
density and pressure of the scalar field and radiation by 



p = p v + p 7 , 



(13) 



where P 7 is the radiation pressure. The energy continuity 
equations for the energy density of the scalar field p^ and 
radiation p 7 in the background are 
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(14) 
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where Q v and Q 7 indicate the background values of the 
transfer coefficient Q v and Q 7 , respectively. 



C« = H 

c 7 1} = 
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(18) 



Notice that the total curvature perturbation £ can be 
expressed as a weighted sum of the single curvature per- 
turbations of the scalar field and radiation fluid as 
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where 
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define the contribution of the scalar field and radiation 
to the total curvature perturbation respectively. On 
the other hand the relative energy density perturbations 
can be characterized in a gauge-invariant manner through 
the difference between the two curvature perturbations 
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The quantity is usually called the entropy or isocur- 
vature perturbation [30]. 



III. GAUGE-INVARIANT PERTURBATIONS AT 
FIRST-ORDER 

The primordial adiabatic density perturbation is as- 
sociated with a perturbation in the spatial curvature ip 
and it is usually characterized in a gauge-invariant man- 
ner by the curvature perturbation on hypersurfaces of 
uniform total density p, usually indicated with £. In- 
deed, both the density perturbations, Sp, Sp^ and 5p 1 , 
and the curvature perturbation, ip, are in general gauge- 
dependent. Specifically they depend upon the chosen 
time-slicing in an inhomogeneous universe. The curva- 
ture perturbation on fixed time hypersurfaces is a gauge- 
dependent quantity: after an arbitrary linear coordinate 
transformation, r — > r + St, it transforms at first-order 
as ip^ — > ip^ + H8t. For a scalar quantity, such as 
the energy density, the corresponding transformation is, 
at first-order, 5^p — ► 5^p — p'Sr. However a gauge- 
invariant combination £ can be constructed which de- 
scribes the density perturbation on uniform curvature 
slices or, equivalently the curvature of uniform density 
slices [7] 

C^=-^-H^. (16) 
P 1 

Similarly it is possible to define the curvature perturba- 
tions Q associated with each individual energy density 
components, which to linear order are defined as [25] 



A. Evolution of first-order curvature perturbations 
on large scales 

The equations of motion for the curvature perturba- 
tions Cp 1 ^ and can be obtained perturbing at first 
order the continuity energy equations (4) for the scalar 
field and radiation energy densities, including the energy 
transfer. Expanding the transfer coefficients Q v and Q 7 
up to first order in the perturbations around the homo- 
geneous background as 

Q ip = Q v + 5 {1) Q ip , (22) 
Q 7 = Q 7 + <5«Q 7 , (23) 

Eqs. (4) give - on wavelengths larger than the horizon 
scale - t 

sv>pf v + m (sw Pv + <5 (1) P V ) - 3 [p v + P v ) v> (1) ' 

= aQ v 4>^+a5^Q v , (24) 
<f<V 7 + m (<5«p 7 + <5 (1) P 7 ) - 3 { Pl + P- ( ) v> (1) ' 
= aQ 7 (1) + a< 5 (1) Q 7 . (25) 



^Here and in the following , as in Refs. [33,35], we neglect 
gradient terms which, upon integration over time, may give 
rise to non-local operators. However, note that these gradient 
terms will not affect the gravitational potential bispectrum on 
large scales. 
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Notice that the oscillating scalar field and radiation have 
fixed equations of state with 5^P V = and S^P y = 
(which correspond to vanishing intrinsic non- 
adiabatic pressure perturbations). Using the perturbed 
(0 — 0)-componcnt of Einstein's equations for super- 
horizon wavelengths ipW + n<p^ = — -j^TT 2 ' we can 
rewrite Eqs. (24) and (25) in terms of the gauge-invariant 



curvature perturbations (jp and [25] 
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2pV P'i P' 



(26) 



(27) 



where 5^Q 7 = —S^Q V from the constraint in Eq (2). 
If the energy transfer coefficients Q v and Q 1 are given 
in terms of the decay rate T as in Eq. (5), the first order 
perturbation are respectively 



5^Q 1 = T5^p ip + 5^T Pv , 



(28) 
(29) 



where notice in particular that we have allowed for a 
perturbation in the decay rate T, 



r(r,x) =r(r) + <5«r(T,x). 



(30) 



Such perturbations in the inflaton decay rate are indeed 
the key feature of the "inhomogeneous reheating" sce- 
nario [17-19]. In fact from now on we shall consider the 
background value T of the decay rate as constant in time, 
r « r» as this is the case for the the standard case of 
inflation and the inhomogeneous reheating mechanism. 
In such a case S^T is automatically gauge- invariant (for 
a gauge-invariant generalization in the case of T' =/= 0, 
see Ref. [29]). Plugging the expressions (28-29) into 
Eqs. (26-27), and using Eq. (19), the first order curva- 
ture perturbations for the scalar field and radiation obey 
on large scales 



C ( l )' = ^^^(c«-c«)+^ (1)r > 

2 p' v P V v > p' v 



(31) 



C (D' 
s 7 



P'i 



Tp' p f (i-g) (c (1) -C«)+W (1) r 



(32) 



From Eq. (19) it is thus possible to find the equation 
of motion for the total curvature perturbation (p) using 
the evolution of the individual curvature perturbations 
in Eqs. (31) and (32) 



c (1) ' = /'(4 1) -c 7 1) ) + /4 1) ' + (i-/K 7 1) ' 

= Hf(l f) (C« - C 7 X) ) = -Hi (C (1) - #>) • (33) 

Notice that during the decay of the scalar field into the 
radiation fluid, p'^ in Eq. (32) may vanish. So it is con- 
venient to close the system of equations by using the two 
equations (31) and (33) for the evolution of ffl and . 



IV. SECOND-ORDER CURVATURE 
PERTURBATIONS 

We now generalize to second-order in the density per- 
turbations the results of the previous section. In par- 
ticular we obtain an equation of motion on large scales 
for the individual second-order curvature perturbations 
which include also the energy transfer between the scalar 
field and the radiation component. 

As it has been shown in Ref. [33] it is possible to define 
the second-order curvature perturbation on uniform total 
density hypersurfaces by the quantity (up to a gradient 
term) 



C (2) = _ v ,(2) _ H 



8™p 



^S^pS_^p U{1)) 



2 



6( r 



P' 



n— -n'- m- 
p 



(34) 



where the curvature perturbation tp has been expanded 
up to second-order as tp = tp^ + ^tp^ and S^p cor- 
responds to the second-order perturbation in the total 
energy density around the homogeneous background p(r) 

p{t,x) = p(t) + Sp(r,x) 

= p(t)+5^p(t^)+ 1 -5^p(t^). (35) 

The quantity C*- 2 -* is gauge-invariant and, as its first-order 
counterpart defined in Eq. (16), it is sourced on super- 
horizon scales by a second-order non-adiabatic pressure 
perturbation [33]. 

In the standard scenario where the generation of cos- 
mological perturbations is induced by fluctuations of a 
single inflaton field (and there is no curvaton) the evo- 
lution of the perturbations is purely adiabatic, and the 
total curvature perturbation is indeed conserved. In 
Ref. [34] the conserved quantity has been used to fol- 
low the evolution on large scales of the primordial non- 
linearity in the cosmological perturbations from a period 
of inflation to the matter dominated era. On the con- 
trary in the curvaton and inhomogeneous reheating sce- 
narios the total curvature perturbation £( 2 ) evolves on 
large scales due to a non-adiabatic pressure. In a sim- 
ilar manner to the linear order it is possible to follow 
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the evolution of C,^ through the evolution of the density 
perturbations in the scalar field and radiation. 

Let us introduce now the curvature perturbations Q 
at second-order for each individual component. Such 
quantities will be given by the same formula as Eq. (34) 
relatively to each energy density pi 



C 



(2) 



(2) 
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Pi 
6^ Pi 



Pi 



H^-H 1 

P l 



:IH 2 ) . (36) 



Since the quantities and Ci^ are gauge-invariant, we 
choose to work in the spatially flat gauge V''- 1 -' = V-^ 2 " 1 = 
if not otherwise specified. Note that from Eqs. (17) and 

(18) and £ 7 ^ are thus given by 
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(37) 
(38) 



Using Eqs. (37) and (38), the energy continuity equa- 
tions at first order (24) and (25) in the spatially flat 
gauge V» (1) = yield 
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We can thus rewrite the total second-order curvature per- 
turbation C (2) in Eq. (34) as 
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where 



1/3 is the radiation equation of state. Using 



Eqs. (41) and (42) to express the perturbation of the total 
energy density 8^ p one obtains the following expression 
for the total curvature perturbation 

C (2) = /C< 2) + (i - /)C 7 2) 

+ /(!-/)(! + /) (C^-Cf) 2 



Q v 
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The (O-O)-component of Einstein equations in the spa- 
tially flat gauge at first-order ipW = 



4>w 



2 p 



2Wp s 



(44) 



and the explicit expressions of the first order perturbed 
coefficients in terms of the decay rate T, Eqs. (28-29), 
yield 



c (2) = /#> + (i - .ncf + /(i - /)(i + /) (#> - c 7 1} ) 

+ ^/(c«-Cf) 2 -2a^)r^ 

aT 

+ 



Pv 



H^\ P 



(45) 



Eq. (45) is the generalization to second-order in the per- 
turbations of the weighted sum in Eq. (19) and extends 
the expression found in Ref. [35] in the particular case of 
the curvaton scenario, under the sudden decay approxi- 
mation, where the energy transfer was neglected. 



where we have used the background equations (14) 
and (15) to find H^r -H' - 2H 2 = -H 2 (6 - /). 

Following the same procedure the individual curvature 
perturbations defined in Eq. (36) are given by 



4 2) = + [2-3(1 + (#>)' 



-2 a 



+ a- 



Q' v a Q v p' 



Hp' 2 Hp' p 



(c»>) 2 , 



(41) 



A. Large-scale evolution of second-order curvature 
perturbations 

In this section we give the equations of motion on large 
scales for the individual second-order curvature pertur- 
bations c}p^ and £ 7 2 \ The energy transfer coefficients Q v 
and Qj in Eqs. (3) perturbed at second-order around the 
homogeneous backgrounds are given by 



Q V = Q V + S (1) Q V + ^Q V , 
Q 7 = Q 7 + ,5«Q 7 ,+i<5( 2 )Q 7 . 



(46) 
(47) 
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Note that from Eq. (2) it follows <5^Q 7 = -S^Qy. Thus 
the energy continuity equations (4) perturbed at second- 
order give on large scales 



6?p 



(1)' 



S^ Pv + 5^P V + 2(/v + P v )ip 
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a 6WQ V + aQ^ - aQ v (V 1 ))' + 2a^5^Q v , (48) 
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a Q 7 p' 



(52) 



Hp' 7 2Hp' 1 pJ\ 

where we have used that w v — and w 7 = 1/3. 

Equations (51) and (52) allow to follow the time- 
evolution of the gauge-invariant curvature perturbation 
at second-order. 



<5< 2 V + m (V 2 )p 7 + <5^p 7 ) - 3(p 7 + p^y 
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5 {1) Pl + 5 {1) P y + 2( Pl + P- f )^ 



a 6 { ®Q 1 + aQ 7 (2) - aQ 7 (<^ (1) ) 2 + 2a#< 1 W 1 >Q 7 , (49) 

where (/>( 2 ) is the second-order perturbation in the gravi- 
tational potential = + ±</> (2) . Note that Eqs. (48) 
and (49) hold in a generic gauge. We can now recast 
such equations in terms of the gauge-invariant curvature 

(2) (2) 

perturbations Q> and Q in a straightforward way by 
choosing the spatially flat gauge tp^ — tp^ = 0. 

The (O-O)-component of Einstein equations in the spa- 
tially flat gauge at first order is given by Eq. (44), and 
at second-order on large scales it reads (see Eqs. (A. 39) 
in Ref. [32]) 



1 5< 2 > , 



+ 4(0«) 2 . 



(50) 



Using Eqs. (44) and (50) with the expressions (41-42) we 
find from the energy continuity equations (48-49) that the 
individual second-order curvature perturbations obey on 
large scales 



c< 2) ' = 

aH " 

V 

iaQ v — 

rip 



S^Q v -^S^pJ+Q v 



p' (5V>p v 8V>p 



'2p 



2a ^(l)Q^(l)_ 2C (l) C (l)' 



p' v 

''W v 2 Hp' v P 



P' v 
a Q v p' 



V. EVOLUTION OF SECOND-ORDER 
COSMOLOGICAL PERTURBATIONS IN THE 
VARIOUS SCENARIOS 

The results contained in the previous section can be 
used in order to study the evolution of the second-order 
curvature perturbations during the reheating phase after 
a period of standard single field inflation, and in the alter- 
native scenarios for the generation of the primordial adi- 
abatic perturbations which have been recently proposed, 
namely the curvaton scenario [13-15] and the inhomoge- 
neous reheating [17-19]. In fact in each of these scenarios 
a scalar field is oscillating around the minimum of its po- 
tential and eventually it decays into radiation. The evo- 
lution at second-order of the curvature perturbations is 
necessary in order to follow the non-linearity of the cos- 
mological perturbations and thus to accurately compute 
the level of non-Gaussianity including all the relevant 
second-order effects. We shall first extend some of the 
results previously obtained in Ref. [35] for the standard 
scenario and the curvaton scenario. Then we will present 
the calculation of the second-order total curvature per- 
turbation produced in the "inhomogeneous reheating" by 
the spatial fluctuations of the inflaton decay rate V. 



A. Standard Scenario 



and 

aH 
— 3aQ 7 



<5< 2 

n 

p' 



5& 



Pi 



P'-r 



2 2a^>)Q 7 ^)-4C 7 1 )C 7 1) ' 
Pi 



<5< 2 



In the standard inflationary scenario the observed den- 
sity perturbations are due to the fluctuations of the in- 
flaton field itself. When inflation ends, the inflaton os- 
cillates about the minimum of its potential and decays, 
thereby reheating the universe. In such a standard sce- 
nario the inflaton decay rate has no spatial fluctuations. 
Here we want to show that in fact during the reheating 

(51) 

phase the curvature perturbations and £( 2 ) do remain 
constant on superhorizon scales. Eq. (33) and Eq. (31) 
with S^T = now read 



(53) 
(54) 



C d)' = _ Hf ( C (D _ <gr 
4 iy = f^7(c (1) -C«)- 



At the beginning of the reheating phase, after the end 
of inflation, the total curvature perturbation is initially 
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given by the curvature perturbation of the inflaton fluc- 
tuations Ci^ = C£L- Therefore Eqs. (53) and (54) show 

that during the reheating phase = = (^J n is a 
fixed-point of the time-evolution. This result has been 
shown at first-order in [29] and extended to second-order 
in the perturbations in Ref. [35] under the sudden decay 
approximation. Under such an approximation the indi- 
vidual energy density perturbations (and hence the cor- 
responding curvature perturbations) are separately con- 
served until the decay of the scalar field, which amounts 
to saying that in the equations for the curvature per- 
turbations Eqs. (51) and (52) one can drop the energy 
transfer triggered by the decay rate aT /H <C 1. Going 
beyond the sudden decay approximation, the first order 
results 4 1} = C (1) in Eq. (45) yield 

C (2) =/4 2) + (l-/)4 2) - (55) 

Deriving this expression and using Eqs. (51) and (52) 
with S^T — and = the equation of motion for 
£( 2 ) on large scales reads 

C (2) ' = -W/(C (2) -4 2) ) • (56) 

In the same way as at first order from Eqs. (55) and (56) 
it follows that the second-order curvature perturbation 
£( 2 ) remains constant on large scales during the reheating 
phase, being given at the end of inflation by the curvature 

(2) (2) 

perturbation in the inflaton field Q n ' = Q j n . 

B. Curvaton scenario 

The evolution of cosmological perturbations in the cur- 
vaton scenario has been studied at second-order in Ref. 
[35] and here we only briefly summarize the main results. 
In the curvaton scenario the final curvature perturbations 
are produced from an initial isocurvature perturbation 
associated to the quantum fluctuations of a light scalar 
field (other than the inflaton), the curvaton a, whose 
energy density is negligible during inflation. The curva- 
ton isocurvature perturbations are transformed into adia- 
batic ones when the curvaton decays into radiation much 
after the end of inflation [14,16]. During inflation, since 
no curvature perturbation is produced = = 0. 
After the end of inflation, the curvaton starts to oscillate 
when its mass is of the order of the Hubble rate and the 
first-order total curvature perturbation can be expressed 
as 

c (1) = /cS 1) + (i - , (57) 

where / = p' a /p' and we have set $p = consistently 
with the curvaton hypothesis. When the curvaton en- 
ergy density is subdominant, the density perturbation in 
the curvaton field gives a negligible contribution to 



the total curvature perturbation, thus corresponding to 
an isocurvature (or entropy) perturbation. On the other 
hand during the oscillations p a oc a -3 increases with re- 
spect to the energy density of radiation p 1 oc a~ 4 , and 
the perturbations in the curvaton field are then converted 
into the curvature perturbation. After the decay of the 
curvaton, during the conventional radiation and matter 
dominated eras the total curvature perturbation will re- 
main constant on superhorizon scales at a value which, 
in the sudden decay approximation, is given by 

c {i) =f D e\ (58) 

where D stands for the epoch of the curvaton decay. Go- 
ing beyond the sudden decay approximation it is possible 
to introduce a transfer parameter r defined as [16,25] 

C (1) = r(W , (59 ) 

where £^ j s evaluated well after the epoch of the cur- 
vaton decay and C,^ is evaluated well before this epoch. 
The numerical study of the coupled perturbation equa- 
tions has been performed in Ref. [25] showing that the 
sudden decay approximation is exact when the curvaton 
dominates the energy density before it decays (r = 1), 
while in the opposite case 

-<j) D <60 » 

At second-order, using Eqs. (51) and (52) under the 
sudden-decay approximation, the individual curvature 
perturbations and ffl are separately conserved and 
the total curvature perturbation (( 2 ) reads [35] 

C (2) = /C (2) + (1 _ /)C (2) 

+ /(l-/)(l + /)(tf ) -# ) ) 2 • (61) 

The second-order curvature perturbation in the standard 
radiation or matter eras will remain constant on super- 
horizon scales and, in the sudden decay approximation, 
it is thus given by the quantity in Eq. (61) evaluated at 
the epoch of the curvaton decay 

C (2) =rC (2) +r(l _ r2) ( C (l)) 2 ; (62) 

where we have used the curvaton hypothesis that the 
curvature perturbation in the radiation produced at the 
end of inflation is negligible so that « and ffl « 0. 

Taking into account that = | (C^) [35], one finally 

obtains the curvature perturbation during the standard 
radiation or matter dominated eras 

C (2) =rg-r 2 )( C «) 2 . (63) 
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C. Inhomogeneous reheating Scenario: 8F =fc 

Recently, another mechanism for the generation of cos- 
mological perturbations has been proposed [17-19]. It 
acts during the reheating stage after inflation and it was 
dubbed the "inhomogeneous reheating" mechanism in 
Ref. [19]. The coupling of the inflaton to normal mat- 
ter may be determined by the vacuum expectation value 
of fields x's of the underlying theory. If those fields are 
light during inflation, fluctuations 5\ ~ where H is 
the Hubble rate during inflation, are left imprinted on 
super-horizon scales. These perturbations lead to spatial 
fluctuations in the decay rate T of the inflaton field to or- 
dinary matter, ^ ~ ^ , causing adiabatic perturbations 
in the final reheating temperature in different regions of 
the universe. Using cosmic time dt — adr the first order 
equation (31) for on large scales now reads 



* 2 p v p \ v J p v 



(64) 



where H = a/ a is the Hubble parameter in cosmic time 
and the dots stand for differcntation with respect to cos- 
mic time. Perturbing the energy transfer Q v = —Tp v up 
to second-order and expanding the decay rate as 

r = r» + 5r = r„ + (5 (1) r + i(5( 2 )r, (65) 

it follows from Eqs. (46) and (65) 

5 {2) Qv = -p v 8V>T - r.SM Pv - 2<5«r<5( V . (66) 
Plugging Eq. (66) into Eq. (51), the equation of motion 

(2) 

on large scales for the curvature perturbation Q, ' allow- 
ing for possible fluctuations of the decay rate 5^T and 
S^T turns out to be 



H 

p'<p 

~it~p 

2H 

p'<p 



Pv 



3r*/v-^ (1)2 

Pv 



+ 



V Pf Pv 

H V P V P 



S {1) Pv 



Pv 



(67) 



where we have used the fact that the decay rate T in the 
scenario under consideration remains constant. 

We shall now solve Eqs. (64) and (67) under a "mixed- 
sudden decay approximation" . We shall treat the pres- 
sureless scalar field fluid and radiation fluid as if they 
are not interacting until the decay of the inflaton when 
r « H. Since at the beginning of the reheating phase 
the energy density in radiation is negligible this means 



that / = p'^/p ~ 1 and there is indeed only a single 

fluid with, from Eq. (19), C (1) « and ffl w 0. In 
fact under such an approximation we can neglect all the 
terms proportional to the decay rate T, but we allow for 
the spatial fluctuations of the decay rate. Thus the first 
order equation (64) reads 



(68) 



where we have used p' v = —3Hp v in the sudden decay 
approximation. Integration over time yields 



Sip 



t 



2 fi')r 

9 H 



C (D 



(69) 



where we have used that during the oscillations of the 
scalar field which dominates the energy density H = ^. 
The inhomogeneous reheating mechanism produces at 
linear level a gravitational potential which after the re- 
heating phase, in the radiation dominated epoch, is given 
by (in the longitudinal gauge) [17] 



v 9 r* 



(70) 



During the radiation dominated epoch the usual relation 
between the Bardeen potential and the curvature pertur- 
bation is tp^ — — §C^\ and thus from Eq. (69) we can 
set the ratio T*/Hd = § at the time of the inflaton decay 
in order to reproduce the result (70) of Ref. [17]. Thus 
from Eq. (69) it follows that the value of is 



C (D ^ _ 



ijWt 
6 r* 



(71) 



We can use this result in order to solve the second-order 
equation (67) , which under the sudden decay approxima- 
tion and using Eq. (44) reads 



C«^-|«Wr-4Wr-2(c«<< 1 )) 



H 



Sip 



(72) 



The fluctuations 5T = S^T + ^5^T indeed depends 
on the underlying physics for the coupling of the inflaton 
field to the other scalar ficld(s) \- Let us take for example 
r(t, x) oc x 2 (i; x )- If the scalar field x is ver Y light, its 
homogeneous value can be treated as constant x{t) ~ X* 
and during inflation quantum fluctuations S^x around 
its homogeneous value \* are left imprinted on super- 
horizon scales. Therefore non-linear fluctuations (S^'x) 
of the decay rate T are produced as well 



rcx X 2 = x* + 2x*<s (1) x+(<5 (1) xy 

From Eqs. (65) and (73) it follows 



(73) 
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p)r = 2 P) x 



5 (2) r /j(i) v \ 2 i /^(i)r N 2 



= 2 



2 v r* 



Using the first order solution ffl = -§<S (1 )r and Eq. (74) 
in Eq. (72), the evolution of (jp on large scales is 

i ? ),__i_( i( .)r) 2 + I( S (.)r)%-2(c W ,) 

Integration over time yields 



r (2) = 



(>r)' 



6 r* 

!A a _2 xl) «Wr 



(c«)-|c 



(76) 



At the time of inflaton decay Y/Hd — f, and since 



4 ' 
.(2) 



H = 2/3 1, it follows t D = 1/2IV Thus Q ' in Eq. (76) 
evaluated at the time to of inflaton decay is 



C ( 2 )._i^V_ c a) 2 _i C (^ (i)r 

V - 24 v r* / ^ 2 V 



(77) 



Using Eq. (71) we finally find that the total curvature 
perturbation £( 2 ) in the sudden decay approximation is 
given by 



c (2) -c< 2) 4(c (1) ) 2 . 



(78) 



VI. SECOND-ORDER TEMPERATURE 
FLUCTUATIONS ON LARGE SCALES 

The goal of this section is to provide the expression 
for the second-order temperature fluctuations on large 
scales which will allow the exact definition of the non- 
linear parameter /nl- From now on, we will adopt the 
Poisson gauge [36] which is defined by the condition 
u> = x = Xi = 0- Then, one scalar degree of freedom 
is eliminated from g oi and one scalar and two vector de- 
grees of freedom from gij. This gauge generalizes the 
so-called longitudinal gauge to include vector and tensor 

(2) 

modes and contains a solenoidal vector w\ . 

The second-order expression for the temperature fluc- 
tuation field in the Poisson gauge has been obtained in 
Ref. [37], by implementing the general formalism intro- 
duced in Ref. [38]. We are interested here in the large- 
scale limit of that expression, which allows to unam- 
biguously define the primordial non-Gaussian contribu- 
tion. Keeping only the large-scale limit of the linear and 



second-order terms in Eqs.(2.27) and (2.28) of Ref. [37], 
we obtain 



AT 
(74) — 



+ v j £ '£ , 
(79) 



where <f>£ 



tfp + \<j> £ 2 ^ is the lapse perturbations at 

emission, T£ = + \r~p is the intrinsic fractional tem- 
perature fluctuation at emission, T£ = AT/T\g. Let us 
recall that, at linear order cj)^ = tp^ . In Eq. (79) we 
dropped all those terms which represent integrated contri- 
butions such as Integrated Sachs- Wolfe, Rees-Sciama and 
more complicated second-order integrated effects [39]. 
Indeed, the non-linear parameter /nl as introduced in 
[27,40] singles out the large-scale part of the second-order 
CMB anisotropies. One should be able to distinguish the 
secondary integrated terms from the large-scale effects 
thanks to their specific angular scale dependence. For the 
very same reason, we disregarded gravitational-lensing, 
time-delay, Doppler terms and all those second-order ef- 
fects which are characterized by a high-£ harmonic con- 
tent. We finally dropped contributions at the observer 
position, which only modify the mono-pole term. 

To obtain the intrinsic anisotropy in the photon tem- 
perature, we can expand the photon energy-density p 1 oc 
T 4 up to second-order and write 



1 S^ Pl 
4 p 7 



3 Pr, 



where p 1 is the mean photon energy-density, and 



-(2) 



1«5< 2 



Pi 



Pi 



3 

16 



p 1 



Pi 



(80) 



(81) 



Next, we need to relate the photon energy-density fluc- 
tuation to the lapse perturbation, which we can easily 
do by extending the standard adiabaticity condition to 
second-order. At first-order the adiabaticity condition 



reads £ 



(i) 



Q and 



we obtain 

±0- {1) Pr, 



Pi 



Pi 



3 Pm 



(82) 



where p m is the average energy-density of the matter 
component. At second-order the adiabaticity condition 

(2) (2) 

imposes Cm — Q ■ From Eq. (36) applied to matter 
and radiation we find 



Pi 



Pi 



4 ^ 

3 Pm 



5™, 



(83) 



In the large-scale limit, the energy constraint, Eqs. (3.3) 
and (4.4) of Ref. [34] in the matter dominated era, yields 



5W, 



-2^ (1) 



(84) 
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and 



-20< 2 > + 8 (V 1 )) 2 . 



We finally obtain the fundamental relation^ 



AT 



(i) 



+ 2 



(85) 



(86) 



From Eq. (86), it is clear that the expression for the 
second-order temperature fluctuations is not a simple 



(i) 



extension of the first-order Sachs- Wolfe effect — - 
\i$P to second-order since it receives a correction pro- 
vided by the term — | ^V^) ■ 

If we express the lapse function at second-order as a 
general convolution (see, e.g., Ref. [32]) 

4> = 4> W + l -^)=^)+f^^)) 2 , (87) 

up to a constant offset, from Eq. (86) we can define 
the true non-linearity parameter /nl which is the quan- 
tity actually measurable by high-resolution CMB exper- 
iments, after properly subtracting instrumental noise, 
foreground contributions and small-scale second-order 
terms. Therefore, we find 



/nl - /] 



NL 



(88) 



We warn the reader that this is the quantity which en- 
ters in the determination of higher-order statistics (as 
the bispectrum of the temperature anisotropies) and to 
which the phenomenological study performed in Ref. [27] 
applies. A number of present and future CMB exper- 
iments, such as WMAP [26] and Planck, have enough 
resolution to either constrain or detect non-Gaussianity 
of CMB anisotropy data parametrized by /nl with high 
precision [27]. 



*One can easily obtain a gauge- invariant definition of the 
temperature fluctuations ^f- in terms of the second-order 
gauge-invariant lapse function (on large scales) 



2 1 V ,W' +2 ^( 1 )) W 1 '" 



„ a nY a la 
+ 5—uj w + — + 




+ 2^' ( 2 0« W X) ') W 2) ' + ^a (2) , 



where 



a< 2 > = ^+3^^' 



+ V" 



-#( 1 >fe< 1 »-J 1 )'9, w ( 1 »-2^ 1 )fe< 1 ) , 



VII. NON-GAUSSIANITY IN THE VARIOUS 
SCENARIOS 

In each of the scenarios described in the previous sec- 
tion it is possible to calculate the level of non-linearity 
in the gravitational potential 4> = <j>^ + \<f>^ (or %p) in 
the longitudinal (Poisson) gauge [36]. For instance the 
gravitational potential 4> can be expressed in momentum 
space as (up to a momentum-dependent function which 
has to be added to account for the condition {<p) = 0) 

0(k) - ^)(k) + ^3 J d 3 ^ d 3 k 2 5® (ki + k 2 - k) 



(27T) 

x/^ L (k 1 ,k 2 )0( 1 )(k 1 )^ 1 )(k 2 ), 



(89) 



where we have defined an effective "momentum- 
dependent" non-linearity parameter . Here the linear 
lapse function i^ 1 ) = tp^ is a Gaussian random field. 
The gravitational potential bispectrum reads 

(0(kiMk 2 )0(k3)> = {2irf <5 (3) (ki + k 2 + k 3 ) 



2 /* L (ki, k 2 ) V^V^kz) + cyclic 



(90) 



where V<j,{k) is the power-spectrum of the gravitational 
potential. We stress again that /^ L defines the non- 
Gaussianity of the gravitational potential, and it does 
not define the non-Gaussianity level of the CMB temper- 
ature fluctuations, see Eq. (86). 

We now give the predictions for the non-Gaussianity 
in the three scenario considered. 



A. Standard scenario 

In the one-single field model of inflation the initially 
tiny non-linearity in the cosmological perturbations gen- 
erated during the inflationary epoch [32,41] gets en- 
hanced in the post-inflationary stages giving rise to a 
well-defined prediction for the non-linearity in the grav- 
itational potentials. In Ref. [34] it has been shown how 
to calculate the second-order gravitational potential at 
second-order from an inflationary epoch to the ra- 
diation and matter dominated epochs by exploiting the 
conservation on large scales of the curvature perturba- 
tion during the inflationary stage and the radia- 
tion/matter phases. 

As we have proved in Section V.A, the curvature per- 
turbation C^ 2 ) is conserved on large scales even during 
the reheating phase after inflation. Using the explicit ex- 
pression for £( 2 ) , the second-order energy constraint (Eqs. 
(4.7) and (4.2) of Ref. [34]) and the traceless part of the 
(i-j)-components of Einstein's equations at second order 
(Eq. (43) of Ref. [35]) we obtain, for perturbations re- 
entering the horizon during the matter-dominated era§, 



Notice that in the Poisson gauge = <^' 2 - ) 



§ Note that this result coincides with the one given in Ref. 
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where 



/* L (k 1) k 2 )~--+ 5 (ki,k 2 ) 



, x ki k 2 / ki • k 2 

ff(ki,k 2 ) = ,„ (1 + 3 



k 2 



k 2 



(91) 



(92) 



where k = ki + k 2 . Notice that in the final bispectrum 
expression, the diverging terms arising from the infrared 
behaviour of /^ L (ki,k 2 ) are automatically regularized 
once the monopole term is subtracted from the defini- 
tion of <j> (by requiring that (</>}=0). Using Eq. (88), we 
conclude that in the standard scenario where cosmologi- 
cal perturbations are generated by the inflaton field, the 
value of non-Gaussianity is provided by 



/NL(ki,k 2 ) 



-- + «/(ki,k 2 ) 



(93) 



B. Curvaton scenario 

In Ref. [35] the level of non-Gaussianity in the gravi- 
tational potential <f> has been calculated following in the 
sudden decay approximation the evolution of the gauge- 
invariant curvature perturbation produced by the 
initial isocurvature perturbations in the curvaton field a. 
In the curvaton scenario we find 



t 
J] 



NL 



7 5 5 
- + -r 

6 6 4r 



+ <?(ki,k 2 ) 



(94) 



Using expression (88), we find that the level of non- 
Gaussianity in the curvaton scenario is provided by 



NL 



1 5 5 
- + -r 

3 6 4r 



+ 0(ki,k 2 ). 



(95) 



C. Inhomogeneous reheating scenario 



C (2) = _^(2) + 



1 6™, 



(2) 



+ 



5 
9 

20 



5^1 



(96) 



where in the last step we have used that on large scales 
8 w p/p = -2-0W in the Poisson gauge [34]. Eq. (96) 
combined with Eq. (78) which gives the constant value 
on large scales of the curvature perturbation during 
the matter dominated era, yields 



v 3 p 6 



(97) 



where we have used the usual relation between the cur- 
vature perturbation and the Bardeen potential C}- 1 ' = 
— during the matter dominated era. From the 

(0-0) and (i-j)-components of Einstein equations, see 
Eqs. (A.39) and (A.42-43) in Ref. [32], the following re- 
lations hold on large scales during the matter dominated 
epoch 



= -|(^ (1) ) + yV" 2 (V 1 W 1) ) 

- 10 V- 4 (d l d 3 • (98) 



We then read the non-linearity parameter for the gravi- 
tational potential <f> — (j)^ + \4>^ 



.^ L = i +. 9 (k 1 ,k 2 ). 



(99) 



Using expression (88), we find that the level of non- 
Gaussianity in the inhomogeneous reheating scenario 
where T cx x 2 is provided by 



/nl = --^ +s(ki,k 2 ) . 



(100) 



Using the technique developed in Refs. [34,35] we can 
now calculate the non linear parameter /^ L in the inho- 
mogeneous reheating scenario. We can switch from the 
spatially flat gauge to the Poisson gauge, since the results 
obtained in the previous section are gauge-invariant, in- 
volving the curvature perturbations. This is evident, for 
example, from Eq. (78). 

During the matter dominated era from Eq. (34) it turns 
out that [34] 



[34] up to an integrated over time gradient term which, how- 
ever, is fully negligible when evaluating the bispectrum of the 
gravitational potential on large scales. 



D. Some comments 

At this point, we would like to make some comments on 
the determination of the non-linear parameter performed 
for the curvaton scenario in Ref. [16] and for the inhomo- 
geneous scenario in Refs. [17,18]. Our results differ from 
those previously obtained for several reasons. First, the 
value of /nl has been obtained in Refs. [16-18] by extend- 
ing at second-order the first-order relation (f>^ — — fC^ • 
This procedure is not correct since this linear relation is 
lost at second-order. Secondly, the value of has been 
obtained for the curvaton scenario in Ref. [16] by merely 
expanding the energy-density fluctuation p at second- 
order and in the inhomogeneous scenario [17,18] by ex- 
panding the decay rate T at second-order. However, 
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the gauge-invariant curvature perturbation contains ad- 
ditional second-order terms of the form (first — order) 2 , 
see Eq. (34). Finally, the value of the non-linear pa- 
rameter /nl has to be defined through the exact relation 
between the temperature anisotropies and the (gauge- 
invariant) gravitational potentials, see Eq. (86). 

VIII. CONCLUSIONS 

In this paper we have presented a general gauge- 
invariant formalism to study the evolution of curvature 
perturbations at second-order. In particular, we have 
addressed the evolution of the total curvature pertur- 
bation in the three scenarios for the generation of the 
cosmological perturbations on large scales: the standard 
scenario where perturbations are induced by the inflaton 
field and the curvaton and the inhomogeneous scenarios 
where the curvature perturbation is produced by initial 
isocurvature perturbations. We have calculated the ex- 
act formula for the second-order temperature fluctuations 
on large scales extending the well-known expression for 
the Sachs- Wolfe effect at first-order. We have provided 
the proper definition of the non-linear parameter enter- 
ing in the determination of the bispectrum of tempera- 
ture anisotropies, thus clarifying what is the quantity to 
be used in the analysis aimed to look for non-Gaussian 
properties in the CMB anisotropies such as the one given 
in Ref. [27]. Finally, we have computed the values of non- 
Gaussianity in the various scenarios for the generation of 
the cosmological perturbations. According to Ref. [27], 
the minimum value of |/nl| that will become detectable 
from the analysis of WMAP and Planck data, after prop- 
erly subtracting detector noise and foreground contami- 
nation, is about 20 and 5, respectively, and 3 for an ideal 
experiment [27]. Our findings imply that a future deter- 
mination of non-Gaussianity at a level much larger than 
./nl ~ 5 would favour the curvaton scenario, while the 
value of non-Gaussianity provided by the standard sce- 
nario is below the minimum value ~ 5 advocated in Ref. 
[27]. However, such a lower bound has been obtained by 
analysing the capability of the bispectrum statistics to 
detect non-Gaussianity of the type discussed in this pa- 
per. Alternative statistical estimators based on the mul- 
tivariate distribution function of the spherical harmon- 
ics of CMB maps [42], or wavelet-based analyses [43], 
together with a systematic application of MonteCarlo- 
simulated non-Gaussian CMB maps [44], may allow to 
achieve sensitivity to lower values of /nl- 
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